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Abstract. We use the photoionisation and dust radiative transfer code MOCASSIN to create 
a model of the dwarf irregular galaxy NGC 4449. The best-matching model reproduces the 
global optical emission line fluxes and the observed spectral energy distribution (SED) spanning 
wavelengths from the UV to sub-mm, and requires the bolometric luminosity of 6.25 x 10 9 Lq 
for the underlying stellar component, Md U st/M ga s of 1/680 and Md us t of 2.2 x 10 6 M©. 
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1. Introduction 

The spectral energy distribution (SED) of a galaxy emerges as a combination of the 
underlying stellar emission, the degree of reprocessing taking place in the interstellar 
medium (ISM) and thermal or non-thermal emission due to dynamical or evolutionary 
processes. As a consequence, galaxies emit across the entire range of the electromagnetic 
spectrum and the distribution of the energy output as a function of wavelength provide 
clues to the galaxy's past and present. Numerical models that reproduce the observed 
global properties of a galaxy over a wide range of wavelengths can be useful in under- 
standing the physical and chemical nature of the underlying stellar and dust components. 
Here we present a MOCASSIN model of the dwarf irregular galaxy NGC 4449. 

MOCASSIN (Ercolano et al. 2003, 2005, 2008) is a fully 3-dimensional photoionisa- 
tion and dust radiative transfer code. Recent applications of the code include modelling 
the planetary nebula NGC 6302 (Wright et al. 2011) and studying dust formation by 
SN 2008S (Wesson et al. 2010). Modelling entire galaxies is technically challenging due 
to their large physical sizes and the high level of complexity of the galactic environments. 
Dwarf galaxies are ideal candidates due to their intrinsically small sizes, relatively simple 
star formation histories and the low degree of processing of their interstellar material. 

NGC 4449 (Fig.0 is a metal-poor magellanic-type irregular galaxy. It is not a typical 
dwarf irregular: at the distance of 4.21 Mpc its H I envelope extends to a radius of 14.2 
kpc (H'6), equal to approximately 6R25 (Swaters et al. 2002), and exhibits a counter- 
rotating core. It was selected for this study for its extensive coverage of photometric 



2 



O. L. Karczewski et al. 




Figure 1. SDSS r-band image of NGC 4449. North is up, east is to the left. The image is 
centred at 12 h 28 m ll!l, +44°05 / 37" (J2000). The bar is 2' in length. 

measurements for wavelengths ranging from the UV to sub-mm. Recently, new FIR 
observations were acquired by the Herschel Space Observatory (Pilbratt et al. 2010) as 
part of the Guaranteed Time Key Project: Dwarf Galaxy Survey (PI: Suzanne Madden). 

The emission line fluxes, used to constrain the models, were derived from global spectra, 
acquired with the drift scanning technique, and are representative of the entire galaxy 
(Kobulnicky et al 1999). 

2. Methods and Results 

Our MOCASSIN simulations use 0.5 million grid cells to establish the observed H I 
profile of Swaters et al (2002) within the inner 3.3 kpc (2'7), enclosing the visible part of 
NGC 4449 and about a half of its total H I mass. The geometry, gas density distribution 
and the observed abundances are fixed inputs. Trial families of simulations are run and a 
range of Md us t/M gas , bolometric stellar luminosities L* and input stellar radiation fields 
are tested. A parameter set providing the best match to the shape of the low-resolution 
SED, the global emission line fluxes and the integrated Ka luminosity (Hunter et al 
1999) is expanded into a new, higher resolution parameter space and the procedure 
is iterated until a satisfactory match with all observables is obtained. Fig. [2] gives a 
schematic summary of this method. 

The stellar radiation field is the most crucial input, as it results from the assumed 



trial synthetic 
radiation field 

trial L* 

trial M dust /M gas 

t 

(expand best 
match parameters 

into a new 
parameter space) 



geometry, n(r) 

4 

MOCASSIN 

>l >l >l 

low-res L(Ha) line 
SED Q(H°) fluxes 

i 

(accept) 



HI, optical, Ha 



(compare) broadband photometry, 
global spectra 



Figure 2. A schematic summary of the modelling method used. The shaded regions group 
together the inputs to MOCASSIN (left) the outputs (middle) and the observational constrains 
(right). 
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star formation history. In this work, two continuous episodes of constant star formation 
activity were assumed and an initial (and subsequently expanded, according to Fig. [2| 
set of 45 templates generated with STARBURST99 (Leitherer et al 1999; version 2011- 
01-17) was tested. A small subset of the resulting SEDs is shown in Fig. [3] 

To break the degeneracy it is helpful to ask if both (i) the shape of the observed SED 
and (ii) the total emergent L(Ra) for a particular star formation scenario can be matched 
by varying M^ ust / M gas and/or L*, whose effect on the resulting SED is relatively simple 
to predict. This condition allows for more than 90 per cent of the initial templates to 
be rejected at an early stage. The final best-matching model is selected by inspecting 
the predicted emission line fluxes. For NGC 4449 the best-matching SED is presented in 
Fig. [4] and the derived physical quantities are given in Table [I] 
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Figure 3. A family of predicted SEDs for four trial star formation scenarios assuming two 
episodes of continuous star formation with two different star formation rates. All input parame- 
ters are identical except for the length of the two star formation episodes (see legend for details) . 
The approximate wavelength ranges affected by the various free parameters are indicated. 
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Figure 4. The predicted SED of the best-matching MOCASSIN model of NGC 4449. PAHs are 
not included in this model. The observed SED is based on the archive data (GALEX in the UV, 
Swift in the UV/optical, SDSS in the optical, 2MASS in the NIR) and the published photometry 
(Engelbracht et al. 2008, Hunter et al. 1986, Bottner et al. 2003; at longer wavelengths). In the 
FIR new Herschel PACS and SPIRE measurements are included (Remy et al, in prep.). 
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parameter MOCASSIN model of NGC 4449 



L* 6.25 ±0.25 x 10 9 L© 



M* 1.2 ±0.1 X 10 9 Mrr 



M gas 1.5 ±0.2 x 10 9 M© 

M dust 2.2 ± 0.2 x 10 6 M© 

M dust /M gas 1/680 (1/850 to 1/540) 

dust composition 100% amorphous carbon 

dust grain sizes 0.005-1 /xm, n oc a -3-5 

SF episodes 6 ± 2 Gyr to 120 ± 40 Myr 120 ± 40 Myr to present 

relative M* 120 1 

<SFR> 0.2 to 0.3 M yr" 1 0.06 to 0.12 M yr" 1 



Table 1. Summary of global parameters derived from the best-matching MOCASSIN model 
of NGC 4449 assuming two continuous episodes of star formation. The global emission line 
fluxes are matched to within 20 per cent or better, with the exception of the sulphur lines [S n] 
AA6717,31, which are overpredicted by a factor of four. 

3. Discussion 

The gas and dust masses derived from our MOCASSIN model of NGC 4449 are con- 
sistent with earlier estimates (Hunter et al. 1999; Engelbracht et al. 2008). Our best- 
matching star formation scenario broadly agrees with the presence of an old (3-5 Gyr) 
and a young population of stars, and continuous star formation over the last 1 Gyr 
(Bothun 1986; Martin & Kennicutt 1997; Annibali et al. 2008). However, the predicted 
present-day star formation rate is only one-fourth of the estimate based on the Ha lumi- 
nosity (Hunter et al. 1999). The discrepancy between these two independent estimates 
highlights the differences in the adopted initial mass function (IMF) as well as the limi- 
tations of the two-episode model. 

The presented technique assumes spherical symmetry for NGC 4449 and a two-episode 
star formation history, both of which are great simplifications. Also, at present PAHs are 
not included in the models and, as a consequence, the 4-20 fim range is not correctly 
predicted. Nevertheless, most of the global parameters are predicted to within the ob- 
servational uncertainties, and therefore this technique may prove useful in studying the 
characteristics of dust and the underlying stellar components in other dwarf galaxies. 
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